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Abstract 

Enhancement of optical Kerr nonlinearity for self-action by electro-magnetically induced trans- 
parency in a four-level atomic system including dephasing between the ground states is studied in 
detail by solving the density matrix equations for the atomic levels. We discern three major con- 
tributions, from energy shifts of the ground states induced by the probe light, to the third-order 
susceptibility in the four-level system. In this four-level system with the frequency-degenerate 
probes, quantum interference amongst the three contributions can, not only enhance the third- 
order susceptibility more effectively than in the three-level system with the same characteristic 
parameters, but also make the ratio between its real and imaginary part controllable. Due to de- 
phasing between the two ground states and constructive quantum interference, the most effective 
enhancement generally occurs at an offset that is determined by the atomic transition frequency 
difference and the coupling Rabi frequency. 

PACS numbers: 42.50.Gy; 32.80.Qk; 42.65.-k 



I. INTRODUCTION 



The weak nonlinear response of even the best materials has been a dominant limitation in 
experimental research on quantum nonlinear optics for many years. A number of theoretical 
proposals, including the creation of a two-photon bound state |j] and few-photon quantum 
solitons , have not yet been experimentally realized due to the lack of large-Kerr-nonlinear 
materials. However, recent research surrounding electro-magnetically induced transparency 
(EIT) ||, which uses atomic coherence to reduce absorption, has opened up a completely 
new route to achieving large optical nonlinearity ||, ^ [7|, 0. An EIT medium gener- 
ally possesses two important features: vanishing resonant absorption and, simultaneously, 
a refractive index curve with a very steep gradient 0. These two features can significantly 
enhance the nonlinear interaction strength in multi-level atomic systems. In addition, the 
latter can also significantly reduce the group velocity of a probe light pulse and therefore 
greatly increase the effective interaction time of the pulse with the medium |7|, |K|. These 
features may therefore enable one to use an EIT medium to achieve nonlinear optical pro- 
cesses at very low light intensities, or even at energies of a few photons per atomic cross 
section |7|, ||. Recently, many EIT-enhanced nonlinear phenomena have been observed in 



experiments, including the Kerr effect P, [11], |T^, [13| . Of particular interest to the present 
work, Wang et al. have measured the Kerr nonlinear coefficient for self-phase modulation us- 
ing a three-level system in Rubidium vapor and demonstrated that the nonlinear coefficient 
is indeed enhanced by EIT |14| . 

In addition to the scheme involving three atomic levels, a four-level system is also a 



candidate for the enhancement of the Kerr nonlinearity for self-phase modulation It is 
not clear at present which, the three-level system or the four-level system, provides the most 
advantages for the enhancement of Kerr nonlinearity for self-action (self-phase modulation, 
two-photon absorption) for very weak light. Or alternatively, because a four- level system in 
some sense contains a three-level subsystem, how does the presence of the fourth level effect 
the enhancement of the nonlinearity? In addition, the four-level system considered in Ref. 
fl5| did not include the dephasing between the two ground states. How does the dephasing, 
which is present in all realistic systems, alter the effective enhancement of the nonlinearity? 
In order to answer these two questions, in this paper we analyze a four-level system with 
dephasing between the ground states. First, in section II we quantitatively compare a four- 



2 



level EIT scheme for self-action with the corresponding three-level scheme. We discern that 
there are three major contributions, due to energy shifts of the ground states induced by 
the probe field, to the third-order susceptibility in the four-level system, and we find that 
quantum interference amongst the three contributions can, not only enhance the third-order 
susceptibility more effectively in the four-level system than in the three-level system, but 
also make the ratio between its real and imaginary part controllable. Next, in section III we 
analyze the detailed behavior of the most effectively enhanced nonlinearity in the four-level 
system. As a result we find that in general, the most effective enhancement does not occur 
exactly at the center of the transparency window due to quantum interference and a finite 
dephasing rate. On the contrary, it occurs at an offset that is determined by the atomic 
transition frequency difference and the coupling Rabi frequency. 



II. THEORETICAL MODEL 



We consider the interaction of two light fields, a coupling field and a degenerate probe 
field, with a gas of atoms, as shown in Fig. [TJ. We assume the atoms can be described by 
the four-level atomic scheme. Consider the case when most atoms are in the ground state 
|1), by applying a strong coupling light between the states |2) and |3), we can dramatically 
reduce the resonant absorption for the weak probe light on the transitions |1) <-» |3) and 
1 2) <-» |4) (see Fig. |l|). Under the rotating wave approximation, this four-level system can be 



described by the following density matrix equations in a frame rotating at frequency l<j p ||16 
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1 

Pll = 731P33 + 2^13/031 - ^3/°13] (1) 

% 

P 2 1 = [«(Ai 3 - A 23 ) - 72l]P21 + -[^23P31 + ^24P41 - ^I 3 P23] (2) 
1 % 

P 3 1 = [iA 13 --(73 1 + 73 2 )]p31 + 2[ fi 13(Pll-P33)+^23P2l] (3) 

1 i 

P41 = [«(A 24 + A13 - A 23 ) - -(740 + 7 42 )]P41 + -[^24P21 - ^1 3 P43] (4) 
P22 = 732P33 + 742P44 + -[^23P32 + ^24P42 - ^ 23 P23 ~ ^ 24 P24] (5) 

1 i 

P 32 = [?A 23 - -(731 + 7 32 )]p 3 2 + ^ [ fi 13Pl2 + ^23(P22 ~ P33) ~ ^ 24 P34] (6) 
1 % 

P42 = [^A 24 - -(740 + 7 42 )]p 4 2 + 2^24(P22 - P44) - ^ 23 P43] (7) 

P33 = -(731 + 732)P33 + 2^13Pl3 + ^23P23 - ^13P31 - ^23P32] (8) 

1 i 

P43 = [z(A 24 - A 23 ) - -(731 + 7 32 + 740 + 742)]P43 + 2 t^24P23 ~ ^13P41 ~ ^23P42] (9) 

1 

P44 = -(740 + 742)P44 + -[^ 24 P24 - ^24P42] (10) 



where p^ is the density matrix element, f2i 3 = fii 3 E p0 /h and f2 24 = p 24 i? p0 /^ are complex 
Rabi frequencies for the probe light field amplitude E p0 , and f2 2 3 = ^zE^/h is the complex 
Rabi frequency for the coupling light with electric field amplitude E c0 , where /113 and p 24 
are electric dipole matrix elements. 721 is the dephasing rate between the ground states |1) 
and 1 2) (This was not included in Ref. The detuning angular frequencies are given by 

A13 = u>p — uJsi, A 23 = ijj c — co> 32 , A 24 = Up — U42, where u 3 i, co> 32 , ^42 are the atomic transition 
frequencies. Additionally, we assume the probe light is very weak, and the coupling light is 
on resonance (A 2 3 = 0). Under this assumption, there is a much larger probability of the 
atoms being in the ground state |1) than in other states, i.e., pn ~ 1. Because the evolution 
of the atoms is very fast and the light is normally approximated as a continuous wave, we can 
consider the atoms to be in steady states. For simplicity of discussion, we assume uniform 
decay rates and uniform electric dipole matrix elements, i.e., 731 = 73 2 = 7 4 o = 742 = 7 and 
P13 — P24 — P- 

Due to the degeneracy of the probe frequency, the polarization induced by the probe field 
is a superposition of the two off-diagonal density matrix elements, 

P = n(p 13 p 3 l + P24P42)e-^* + C.C. (11) 
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where n is the atom density. It is worth noting that this linear superposition gives rise 
to quantum interference between the two off-diagonal density matrix elements, because 
the two quantities are complex functions and there is coherence among the atomic states. 
For example, if the transition amplitudes p^i and p^2 are in phase, an effective two-photon 
transition will be enhanced; on the other hand, if the two amplitudes are out of phase, then a 
photon emitted on one transition will be absorbed on the other and the effective two-photon 
transition will be suppressed. It is evident that there is no interference between the two 
off-diagonal density matrix elements in the four-level system with frequency-non-degenerate 
probes ||. Thus this superposition of the two off-diagonal density matrix elements is a 
unique feature of the four-level system with one frequency-degenerate probe. Because the 
probe field is monochromatic, the corresponding susceptibility has the simple form, 



s {x W (E c0 )E p0 + x (NL) (E c0 , E p0 )E 3 p0 ) = 2n(/i 13 p 3 i + /W4 2 ) (12) 

where x is the linear susceptibility and x^ N the nonlinear susceptibility. When E p0 — > 0, 
t^( jvx ) corresponds to the third-order susceptibility x^- The real part of x^ is proportional 
to the Kerr refractive index and the imaginary part of x^ to the two-photon absorption 
coefficient. Here, x does not exist because of the symmetry of the atomic medium. Al- 
though the occupation probabilities of atoms on both the states |2) and |4) are very small, 
the contribution of p^i to the third-order susceptibility is not small, as we will see in next 
section, therefore neglecting the contribution of p^2 to the third-order susceptibility is in- 



correct | i~5fl . To analyze the enhanced optical nonlinearity of the EIT medium, we use a 
numerical method, Gaussian elimination, to solve Eqs. ([?])- (|TI]) in the steady state and then 
extract the first-order and the third-order susceptibilities via the relation (|T2j). 



III. NUMERICAL COMPARISON OF FOUR-LEVEL AND THREE-LEVEL SYS- 
TEMS 

In this section, we numerically compare a four-level EIT system with a three-level EIT 
system. The three-level system we consider is in fact a special case of the present four- 
level one. When the atomic transition frequency difference U31 — CJ42 is very large, if the 
probe frequency u p is close to the atomic transition frequency u 3 x, it will be far detuned 
from the other atomic transition frequency 0042- In this case, any effect from the state |4) is 
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negligible, and the four-level system can be approximated as a three-level system with similar 
parameters. As the three-level system is actually a subsystem of the four-level system, the 
comparison with the three-level system is not only to show the advantage of the four-level 
system in enhancing Kerr nonlinearity, but also to show the contribution of the three-level 
subsystem, i.e., the states |1), |2) and |3), to the third-order susceptibility in the four-level 
system. 

We do the comparison first in terms of third-order susceptibility x ■> an d then in terms 
of the ratio between the third-order susceptibility and the first-order susceptibility, i.e., 
X^ 3 ' /Imx^ ■ In the first comparison, we show that there are three major contributions of light 
shifts of the ground states to x^ i n the four-level system: one of them is actually produced 
within the three-level subsystem, and other two are specific to the four-level system. In 
particular, we show that there is quantum interference, which does not exist in the four- 
level system with frequency- non-degenerate probes ||, amongst the three contributions due 
to atomic coherence. As a result, one can achieve not only larger third-order susceptibility in 
the four-level system than in the three-level system, but also a control of the ratio between 
its real and imaginary part. Then in the final part of the section, we consider the finite 
dephasing rate between the ground states, and further show the advantage of the four- 
level system in effectively enhancing the third-order nonlinear susceptibility with this new 
criterion. 



A. Analysis in terms of x^ 

First, following the idea described in the beginning of this section, we let U31 — U42 = 
— 10 5 7 and the four-level system becomes effectively a three-level system. Any contribution 
of the state |4) to the susceptibility is negligible, so the third-order susceptibility x^ 3 \ as 
shown by the dashed lines in Fig. (a) and (b), only originates from the three states |1), 
1 2) and |3), through the off-diagonal density matrix element p 31 . Its real part has already 
been demonstrated to be much larger than that of the resonantly enhanced susceptibility 
without the coupling field at a finite detuning fll4j| . Here we can understand the unusually 
large magnitude of Rex® in a similar way to Ref. ||: the presence of the weak probe field 
between the states |1) and |3) causes an energy shift of the state |2), which results in an 
effective shift of the linear susceptibility x^ curve (see Fig. ||] (c)). However, here there is 
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an important difference from Ref. ||: these shifts of the curves are not horizontal shifts 
along the frequency axis, but changes of the gradients of the curves. For Rex 1 - 1 ' , the gradient 
variation is equivalent to a small rotation of the solid curve around the zero detuning. Then 
we can understand that although Rex® in the situation of EIT is much larger than the 
resonantly enhanced Rex*- 3 " 1 without the coupling field at a finite detuning, the magnitudes 
of Rex® and Imx® also become zero at the zero detuning, as shown in Fig. [2] (a) and 
(b). Thus the light shift within the three-level subsystem produces a contribution to the 
third-order susceptibility in the four-level system, and the contribution becomes important 
when the probe detuning is finite. 

Next, let us consider the solo contribution of p^i to x at a small atomic transition 
frequency difference u 3 i — U42 plotted by the thin solid lines and the thin dot-dashed lines 
in Fig. m (a) and (b). We can find that although both the thin solid curves and the thin 
dot-dashed curves keep some resemblance of the dashed curves of the three-level system, the 
magnitudes of some parts of the curves are significantly larger than those of the three-level 
system at the same frequency. These variations indicate that the presence of the state |4) has 
important influence on x through p%i. We can understand this effect in the same way as 
Ref. 0, |J: the presence of the weak probe field between the states |2) and |4) leads to another 
energy shift of the state |2), which results in another effective shift of the linear susceptibility 
X®, i.e., the curves experience additional shifts besides the gradient shifts produced within 
the three-level subsystem. For the real part of the susceptibility, the shift is mainly along 
the frequency axis when the atomic transition frequency difference U31 — uo^ 7^ 0. Thus the 
shift results in a finite value of Rex® at the center of the transparency window. 

Next, following Eq. (|i~2"D , we include the contribution of p 42 in the calculation of x®, 
as shown by the thick solid curves and thick dot-dashed curves in Fig. ^| (a) and (b). In 
comparison with the corresponding thin solid curves and thin dot-dashed curves, we can 
find that both Rex'- 3 '' and Imx® are drastically changed once again. In particular, Imx® 
is twice its previous value at zero detuning when oj^i = U42, as shown by the two solid lines 
in Fig. ^ (b). This indicates that the contribution of P42 to x® is not small, and neglecting 
its contribution as in Ref. [HJ does not give us correct values of x®- We can see that the 
shift results in another increase of the magnitude of Rex® at the center of the transparency 
window when the atomic transition frequency difference U31 — CJ42 7^ 0. Similarly, the energy 
shift of the state |2) due to the interaction of the weak probe field with the states |1) and 
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1 3) leads to another effective shift of x^ through p 42 . 

So far we have found three major contributions of light shifts to the third-order suscep- 
tibility in the four-level EIT system. However, the most important feature is that there 
is quantum interference amongst the three contributions because of the existence of the 
atomic coherence among the atomic levels. The linear superposition of p^i and P42 in Eq. 
(|12|) shows explicitly that the interference can occur. Through comparison of the curves of 
-^-(3) f or different contributions in Fig. ^| (a) and (b), we can also discern the variation of the 
susceptibility due to the interference: the magnitude of x^ ma y be much increased at cer- 
tain detuning due to constructive interference, but much reduced at another detuning due to 
destructive interference. For example, for A 13 »s 0.257, the magnitude of Rex^, plotted by 
the thick dot-dashed curve in Fig. (a), is increased in comparison with the thin dot-dashed 
curve and the dashed curve at the same frequency due to constructive interference; however, 
the magnitude of Imx^ vanishes at the same detuning due to destructive interference, as 
shown in Fig. |2] (b). In the four- level system with frequency-non-degenerate probes, there 
is no such quantum interference effect ||. This quantum interference process is a unique 
feature of the four-level system with degenerate probes that we study in this paper. By 
using this quantum interference we can, not only enhance x^ m the four-level system more 
than that in the three- level system, but also control the ratio between Rex^ and Irax^, 
as shown in Fig. |3|. Comparing the thin curves and the corresponding thick curves, we find 
that a small probe detuning can produce some zero points for Imx^ and change the ratio 
dramatically. If, for example, the four levels are magnetic sub-levels in experiments, one 
can use a magnetic field to control the atomic transition frequency difference u^ 31 — uo i2 and 
simultaneously set the probe detuning A31 to achieve expected ratios. Undoubtedly, such 
control of the nonlinear susceptibility will be very useful in a practical design of various 
optical devices. 

B. Analysis in terms of x /IniX 

Because the magnitude of a nonlinear effect per unit length becomes very small at very low 
light intensities, such as when there are few photons per atomic cross section, propagating 
light beams for a long distance in a nonlinear medium is usually a good way to magnify 
the nonlinear effect. However, the distance that light beams can propagate is limited by 
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the linear absorption of the medium, which cannot be zero even under the conditions of 
EIT because normally 721 is not zero. Thus, the ratio between the nonlinear coefficient 
and the linear absorption of the medium is the real criterion for the effectiveness of the 
enhancement of the nonlinear coefficient. We therefore calculate the ratio between the 
third-order susceptibility x^ an d the imaginary part of the first-order susceptibility x^'- 

A = x (3) /Imx {1) (13) 

The results for the four-level system and the three-level system are shown in Figs. H(a), 
(b), (c), and (d). In order to determine the positions of the peaks of A relative to the 
transparency window, we also depict Im%W in Fig. f|(e). 

We find that for the four-level system with a small atomic transition frequency difference 
co> 31 — CU42 = —7, the position of the largest peak (along the detuning axis) of ReA (or ImA) 
is always close to the center of the transparency window independent of the value of the 
Rabi frequency ^23, as shown in Figs. f|(a) and (b). In contrast, the three- level system can 
only produce very small peaks (along the detuning axis) of ReA at the central part of the 
transparency window, as shown in Figs. |](c) and (d). The largest peaks (along the detuning 
axis) of ReA (or —ImA) of the three-level system are not in the central region but very close 
to the two edges of the transparency window, and move away from the center when the 
transparency window becomes wider, as indicated by the arrows in Figs. [|(c) and (d) - 
Wang et al. have already demonstrated this phenomenon in their recent experiments |T7] . 



Because the magnitude of Imp^ 1 ) at the edges is much larger than at the center [see Fig. 
H (e)], the magnitudes of the largest peaks of A of the three-level system are much smaller 
than those of the central peaks of A of the four-level system as seen in Figs. |](a) and (b). In 
this calculation, we assume the dephasing rate 721 = O.OI7. In this case the largest peak for 
the three-level system is more than one order smaller than the largest peak for the four-level 
system. 

From the above comparison in terms of A, we further confirm that the third-order sus- 
ceptibility x is indeed more effectively enhanced in the four-level system, and therefore 
the four-level system has an advantage for the realization of many quantum nonlinear optics 
phenomena. However, this conclusion does not mean the contribution of the three-level 
subsystem to the effectively enhanced A is negligible. In the next section, we will show the 
important influence of the quantum interference amongst the three contributions of light 
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shifts to x 011 the largest peak of A at a finite probe detuning when the transparency 
window becomes wider. 



IV. BEHAVIOR OF THE LARGEST PEAK OF A OF THE FOUR-LEVEL SYS- 
TEM 

In this section, we consider the detailed behavior of the largest peak of A of the four-level 
system with a small atomic transition frequency difference U31 — uj^ and a finite value of the 



dephasing rate 721. Because the dephasing rate 721 is always finite in reality JTy, [Hfl , the 
absorption described by Imx^ 1 ^ is not zero when A13 = O7, and A is not divergent at this 
detuning. Thus, it is important in an implementation of the four-level scheme to determine 
under what condition the most effectively enhanced nonlinear susceptibility occurs under a 
finite dephasing rate 721. In Fig. [| (b), we can easily discern that the value of the largest 
peak (along the detuning axis) of ImA increases when the coupling Rabi frequency 2 3 
increases. By checking more carefully, we find that not only the value of the peak of A but 
also the position of the peak depends on the coupling Rabi frequency. This is depicted in 
Fig. |. 

Fig. |5](a) shows the dependence of the value of the largest peak (along the detuning axis) 
of ImA, i.e., (ImA) pea fc, on the Rabi frequency, S7 2 3, and the atomic transition frequency 
difference u^i — oj^- In this calculation, we assume the dephasing rate 721 = O.OI7. We 
find that as the coupling Rabi frequency, ^3, increases, (ImA) pea fc (for a fixed value of 
^31 — ^42) has a relatively small value at first, increases very rapidly, and then saturates to 
a relatively large value. For a fixed small value of ^23/7, (ImA) pea fc decreases monotonically 
as (co>3i — uj i2 )M decreases from to —3. However, the saturation values for different values 
of uj 3 i — U42 always appear to be the same, and they are always larger than (ImA) pea fc at 
small ^23/7- Although the value of the largest peak of ReA, (ReA) pea fc, is not a monotonic 
function of (0^31 — oJ^/l a t small ^23/7, its behavior is quite similar to that of (ImA) pea fc: 
as the coupling Rabi frequency ^23 increases, (ReA)p ea fc always saturates at the same value, 
which is larger than its value at small ^23/7 in most shown in Fig. |5|(c). The 

behavior of the largest peaks of ReA and ImA indicate that if the most effectively enhanced 
nonlinear susceptibility is desired under a finite 721, then the coupling Rabi frequency 2 3 
should be set as large as possible. This is especially true when the magnitude of the atomic 
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transition frequency difference is not small. 

Next, from Fig. |5](b) [or (d)], we find that for a fixed value of (0031 — uj^/j, when 
(Im\) peak (or (ReA) pea fc) increases and saturates, the detuning (A 13 ) peafc , at which the largest 
peak occurs, also shifts asymptotically from (A 13 ) pea fc = to another finite value. The 
asymptotic value of the detuning is always equal to — (uj 3 i — u^)/?, as shown in Fig. |§(b). 
This value exactly corresponds to the frequency of the probe light being resonant with the 
two-photon transitions between states |1) <-> |4), i.e., 2u p = u^i + uj^- The asymptotic value 
of the detuning at which the largest peak of ReA occurs is always a half line width, 7/2, away 
from that of the largest peak of ImA for the same (co>3i — ui^/l, as shown in Fig. [|(d). This 
behavior of the largest peaks of ReA and ImA indicates that while the transparency window 
becomes wider, the constructive quantum interference amongst the three contributions of 
light shifts becomes more important, and the peak of A therefore shifts to a nonzero probe 
detuning. Thus if a most effectively enhanced nonlinear susceptibility is desired under a 
finite 721, the detuning of the probe light should be set close to the asymptotic values found 
in the above analysis. This is especially true when the magnitude of the atomic transition 
frequency difference is not very small. Even if a certain ratio between Rex^ and Inr^ 3 -' is 
required, the asymptotic values shown above will be important for determining the optimum 
detuning setting. Setting the detuning to zero, i.e., the center of the transparency window, 
might not be the best choice. 

Additionally, an interesting case occurs when CJ31 — U42 = ~1 (or O7). In this case, 
(ReA)p ea fc (or (ImA) pea fc) reaches its saturation value for very small values of ^23- This 
means that the enhancement of the nonlinear susceptibility can occur at very low coupling 
light intensity. 

In the calculations in Fig. ^ and Fig. |^, we assumed the dephasing rate of the ground 
states is 721 = O.OI7. More generally, we can calculate the dependence of the saturation 
values of the largest peaks of A on the dephasing rate 721, as shown in Fig. |j. If the dephasing 
rate 721 is much smaller than the decay rate 7 of the upper levels, then the four-level system 
can enhance the nonlinear susceptibility more effectively, as indicated by A. This means 
that in order to employ the four-level system to more effectively enhance the nonlinearity, 
one needs to reduce the dephasing rate of the ground states as much as possible. 

The analysis of the behavior of the largest peak of A also demonstrates the importance 
of the quantum interference amongst the three contributions of light shifts to when the 
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transparency window is very wide. To implement the four-level system with a finite 721 
to more effectively enhance the nonlinear susceptibility, the coupling Rabi frequency Q23 
should be set as large as possible and simultaneously the detuning of the probe light should 
be set as close as possible to the asymptotic values found in the above analysis. This is 
particularly important when the magnitude of the atomic transition frequency difference is 
not so small. This conclusion is very different from the proposal in Ref. [T5|] due to the fact 
that in that work only the very special case when the dephasing rate vanishes, 721 = 0, was 
considered. 



V. CONCLUSION 



We have studied in detail the third-order susceptibility for self-action of a four-level sys- 
tem, including the dephasing between the two ground states, under the condition of EIT 
by numerically solving the steady-state equations for the atomic density matrix. Through 
comparison of the four-level system with a three-level system with the same characteristic pa- 
rameters, we discerned three major contributions from light shifts to the third-order suscep- 
tibility in the four-level system. In particular, we found that quantum interference amongst 
the three contributions, which does not exist in the four-level system with frequency-non- 
degenerate probes ||, can not only enhance the third-order susceptibility more effectively 
in the four-level system than in the three-level system, but also make the ratio between its 
real part and imaginary part controllable. This unique feature means the four- level system 
has certain advantages for the realization of many quantum nonlinear optics phenomena. In 
implementing this scheme, it is important to note that in general the most effective enhance- 
ment of the nonlinear susceptibility does not occur exactly at the center of the transparency 
window. Instead, due to the constructive quantum interference, and a finite dephasing rate 
between the two ground states, the most effective enhancement occurs at an offset that is 
determined by the atomic transition frequency difference and the coupling Rabi frequency. 
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FIG. 1: Energy levels and optical couplings of the four-level atomic system. to c is the angular 
frequency of the coupling light, and lo p the angular frequency of the degenerate probe light. Direct 
electric-dipole transition between two ground states, |1) and |2), is assumed to be forbidden. 731, 
732, 742 are decay rates from excited states to the ground states. 740 is the decay rate of the state 
1 4) to states other than these four states. A13 and A 24 are detuning frequencies of the degenerate 
probe light. 
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FIG. 2: Probe susceptibilities vs the probe detuning frequency, (a) and (b) show the real and 
imaginary parts of the third-order susceptibility x*- 3 -* for the four-level system in three cases: W31 — 
U42 = — 10 5 7 (dashed lines, an approximation to a three level system), W31 = LO42 (thick solid lines) 
and W31 — ll>42 = 27 (thick dot-dashed lines), %( 3 ) is i n units of 2n / u 4 /(eo^ 3 7 3 )- I n addition, they 
also show the solo contribution of p^i to m the second case (thin solid lines) and the third 
case (thin dot-dashed lines), (c) shows the transparency window: the real part (solid line) and 
imaginary part (dashed line) of the first-order susceptibility x^ both become zero at A13 = O7. 
%W is in units of 2rifi 2 /(eoh'y). In the calculation, 721 = O7 and ^23 = 27. 
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FIG. 3: The third-order susceptibility vs the atomic transition frequency difference W31 — W42 
in the four- level system for two kinds of probe detunings A 13. The real and imaginary parts of 
in units of 2n / u 4 /(eo^ 3 7 3 ) are represented by the solid line and the dashed line, respectively. 
The inset shows the ratio between the real part and the imaginary part of as a function of 
^31 — ^42- In the calculation, A 24 = W31 — ^42, ^23 = 27 and A13 = O7 (thin lines) or O.47 (thick 
lines). 
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FIG. 4: Ratio A vs the probe detuning A13 and the Rabi frequence ^23. (a) and (b) show the 
ratio for the four-level system with U31 — 0042 = — 7- (c) and (d) show the ratio for the three-level 
system, i.e., the approximation of the four- level system with ojsi — U42 = — I0 5 j. (e) shows the 
imaginary part of the first-order susceptibility ■ The ratio A is in units of ^/(fry) 2 , and the 
linear susceptibility is in units of 2n / u 2 /(eo^7)- I n this calculation, 721 = O.OI7. 
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FIG. 5: Value and detuning of the largest peak of A vs the Rabi frequency Q23 and the atomic 
transition frequency difference W31 — U42. (a) and (b) are for the imaginary part of A. (c) and 
(d) are for the real part of A. The curves for (0*31 — W42V7 = 0, —0.5,-1,-2,-3 are plotted by 
dot-dashed lines, dashed lines, dotted lines, thin solid lines and thick solid lines, respectively. The 
ratio A is in units of // 2 /(fry) 2 . In this calculation, the dephasing rate 721 = O.OI7. 
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FIG. 6: Saturation value of the largest peak of A vs the dephasing rate 721. The ratio A is in units 
of /^/(hj) 2 . The solid line is for the real part of A, and the dashed line for the imaginary part of 
A. 
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